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We have studied the effects of electron correlation on Van Vleck susceptibility (xvv) in tran- 
sition metal compounds. A typical crossover behavior is found for the correlation effect on xvv 
as sweeping spin-orbit interaction, A. For a small A, orbital fluctuation plays a dominant role in 
the correlation enhancement of xvv; however, the enhancement rate is rather small. In contrast, 
for an intermediate A, xvv shows a substantial increase, accompanied by the development of 
spin fluctuation. We will discuss the behavior of xvv in association with the results of Knight- 
shift experiments on S^RuCU and an anomalously large magnetic susceptibility observed for 
5d Ir compounds. 
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Multi-orbital systems often show intriguing phenom- 
ena, such as unconventional superconductivity 1, 2 ' and 
colossal magneto-resistance. 3 ' In most cases, orbitals do 
not just complicate a system by increasing the number 
of local degrees of freedom, but they lead to qualitatively 
new features, which are absent in single-orbital systems. 
Among the characteristic properties inherent in multi- 
orbital systems, the role of orbital angular momentum 
and the spin-orbit coupling deserves special attention, 
since they underlie many anomalous phenomena in tran- 
sition metal and rare-earth materials. 

Above all, these two factors give rise to "residual" 
paramagnetic susceptibility, called Van Vleck suscepti- 
bility (xvv)- 4,5 '' It is well known that the Pauli suscep- 
tibility (xp) is proportional to the density of states, and 
it vanishes, as soon as an energy gap opens at the Fermi 
level. Meanwhile, xvv has a residual nature, in the sense 
that it remains finite at zero temperature, even in the 
presence of an energy gap. 

To be more specific, the residual nature of xvv can 
be attributed to the fact that the orbital angular mo- 
mentum or spin-orbit interaction makes magnetization 
a non-conserved quantity. To illustrate this, suppose a 
system in which the magnetization is conserved. In this 
case, the Hamiltonian and magnetization operators are 
commutative, and thus a fixed magnetization quantum 
number can be assigned to each eigenstate of the Hamil- 
tonian. Consequently, for the magnetic susceptibility to 
be finite at zero temperature, the ground state has to 
be replaced with one of the excited states with an in- 
finitesimal magnetic field, which is obviously impossible 
in a gapped system. Meanwhile, for a system where the 
magnetization is not conserved, an infinitesimal magnetic 
field leads to magnetization by mixing the ground-state 
and excited-state wave functions, resulting in a finite sus- 
ceptibility even in the presence of an energy gap. 

In fact, the residual nature of Xvv causes much confu- 
sion in gapped systems. A typical example can be found 
in a Knight-shift experiment to determine the parity of 



a superconducting order parameter. Below the supercon- 
ducting transition temperature T c , only xp decreases, 
responding to spin gap formation, while xvv remains in- 
variant. This means that the invariant Knight shift at 
T c does not necessarily serve as evidence of spin-triplet 
superconductivity, but suggests another possibility that 
a very large xvv masks the decrease in xp to experi- 
mental accuracy. The most famous material that suffers 
from this difficulty may be UPt3, for which a decrease 
in Knight shift has been observed; however, the decrease 
is only 1% of the total Knight shift. 6 ' Similar difficul- 
ties have been reported for conventional vanadium-based 
superconductors. 7 ' Recent Knight-shift experiments on 
Si'2Ru04 also led to the same question. Murakawa and 
coworkers found that the magnetic susceptibility shows 
no change at T c , irrespective of the magnetic field direc- 
tion. 8,9 ' In order to provide a basis for discussing these 
interesting superconductors, it is desirable to clarify the 
Van Vleck susceptibility of multi-orbital systems. 

5c? Ir compounds provide another example in which 
Van Vleck susceptibility plays a crucial role. One of the 
confusing properties common in Ir compounds is their 
unusually large magnetic susceptibility and Wilson ra- 
tio. For example, for EU2R2O7, the magnetic suscepti- 
bility amounts to x ~ 1.0 X 10 _2 emu/mol-Ir, in contrast 
to the rather small specific heat coefficient 7 ~ 8.0mJ/K 2 
mol-Ir, leading to an anomalously large Wilson ratio, 
i?w ~ 90. 10 ' A similar huge paramagnetic susceptibil- 
ity has also been reported for the so-called "hyperkagome 
material" Na^rsOs, where a large residual magnetic sus- 
ceptibility, x ~ 1-0 x 10 _3 emu/mol-Ir was observed, 11 ' 
despite that this compound is an insulator. To recon- 
cile a large paramagnetic susceptibility with a small 7, 
substantial increase in xvv is necessary, since a large 
Xp requires the existence of rich gapless spin excitations, 
which should also contribute to 7. 

Actually, in rare-earth systems, correlation effects on 
Van Vleck susceptibility have been studied by several 
groups. 12 ~ 17 ' Among them, Kontani and Yamada pointed 
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out the general tendency that the correlation enhance- 
ment of xvv is comparable to that of xp. 14 ^ However, it 
has also been reported that the enhancement rate consid- 
erably depends on individual properties of a system, such 
as orbital degeneracy. 15, 17 -* Therefore, it is highly non- 
trivial how electron correlation affects xvv hi transition 
metal compounds, which have quite different characters 
from rare-earth materials. 

In this research, we will study the correlation effect on 
Xvv m transition metal compounds. For this purpose, 
we adopt Si'2Ru04 as a model material, since its simple 
orbital structure is appropriate for establishing a general 
theory, and the xvv of Si'2Ru04 is interesting of its own 
right. Although the spin-orbit coupling of Sr2Ru04 is 
rather small, we also investigate the case with a large 
spin-orbit coupling, in order to gain insights into large 
Xvv in 5d transition metal compounds. Hereafter, we 
set h=k,B = fJ-B = 1 • 

We start with a multi-orbital Hubbard model, which 
takes account of the three ti g orbitals of Sr 2 Ru04 , 

H = H Q + Ht, (1) 
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Here, s = +1(— 1) denotes the up- (down-) spin. The 
indices a, b = 1,2, and 3 represent d yz , d zx , and d xy 
orbitals, respectively. Here, we introduce the abbrevia- 
tion i] = (s, a). We define \yz) = -^=(| + 1) + | - 1)), 

\zx) = -^(| + 1) - | - 1)), \xy) = -75(1+2)-!- 2)), 
where |m) is the eigenstate of the orbital angular mo- 
mentum, with l z \m) = m\m). This convention makes Hq 
real, at the cost of the sign factor in pair hopping terms: 
C, a b = —1, if o, = 2 or b = 2, and ( a b = 1, otherwise. We 
determine the kinetic energy terms in Hq, by the two- 
dimensional tight-binding model: ei(k) = — 2^ cos — 
2t xy cos k y — fi, 62 (k) = —2t xy cos k x — 2t' xy cos k y — \x and 
63(k) = — 2i z (cos k x + cos k y ) — At' z cos k x cos k y — [i. We 
set t z = 1 as a unit of energy, and t xy = 1.5, t' xy = 0.2, 
t', = 0.4 to reproduce the Fermi surface obtained in the 
de Haas- van Alphen experiment. 18 -' The chemical poten- 
tial fjj is controlled so that the system is at 2/3 filling. 

In order to separate xp an d Xvv> we follow the dis- 
cussion by Kontani and Yamada. 14 ) We fix the magnetic 
field parallel to z; then, we define 



Xvv = lim lim x q *M, 



(4) 



X = Xp + Xvv = lim, lim x q z (w), 

q— >0uj— >0 



Xp = X - Xvv, 



(6) 



where the magnetic correlation function Xq Q (w) can be 
obtained from its Matsubara-frequency representation 



Xq=n(3^T)- x ^ (i7rr) 
2 



X%S(iu q ) = jjjdr e-^(M« q (r)M«') (a,a' = x,y,z) 


as Xq"'(w) = lim |8 _ ) . 00 Xq Q '(«^ g -> oj + iS), where M™ = 

Ek, W '(^l^ a l'7 , )4+ q ,^k, I7 ', with M a = P + 2s«, and N 
is the number of sites. The definitions eqs. (4)-(6) guar- 
antee the properties required for Pauli and Van Vleck 
susceptibilities: xp vanishes if an infinitesimal energy gap 
opens at the Fermi level. These definitions are meaningful 
only at zero temperature. However, for numerical calcu- 
lation, we introduce a small temperature, T, and approx- 
imate xvv as xvv ^ X q =o( i7rT ) 
We typically set T = 0.01. 

In order to obtain XqC^g)) we use the dynamical 
mean-field theory (DMFT), 19 ) in which the self-energy 
and the irreducible vertex function are approximated by 
the local ones. In our model, the local Green's function 
can be expressed as a 6 x 6 matrix, G m > (*e p ), which satis- 
fies the self-consistent equations involving the Weiss field, 
Q-qri 1 (*£p)- 19 ^ Here, we obtain the self-energy Y, m i(ie p ) us- 
ing the iterative perturbation theory (IPT): we expand 
the self-energy to the second order of electron interac- 
tions, U, U', and J. This method is valid only for a weak- 
coupling region, however, computationally inexpensive 
and appropriate for studying a wide parameter range. 
The self-energy can be written as 



^mv2 ( r ) — ^mmmvidmrn ( r — ) 



where we have introduced the bare vertex function 

^viV2V3ni = 2(/t;it/2»)3';4 — Iv2ViV3"<h) ■ Then the mag- 
netic correlation function is written as 
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where the two-body Green's function Gq_ mri2ri3rill can be 
obtained from the Bethe-Salpeter equation, 19 ) in which 
we calculate the local irreducible vertex function up to 
the second order of electron interactions as 
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Below, we show our results. Firstly, we will discuss 
Xp and xvv f° r a non-interacting case and A = 0. In 
this case, the definitions of xp an d xvv, e qs. (4)-(6) 
can be naturally extended to a finite temperature. We 
add the Zeeman term Hz = -h J2i, m < (v\ MZ W) c lr, c iv' 
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to the Hamiltonian cq. (1) and obtain xp = 
^E k Xp(k) = ^E k ,a(%^) 2 (-<§ /Mk) -M)), and 
Xvv = ^EkXvv(k) = -iEmIt)^) - ^' 
with a Fermi distribution function, /(x) = e ^ x 1 +1 ■ We 
plot the temperature dependences of xp and xvv m Fig. 
1 (a), xp shows a moderate increase with decreasing tem- 
perature, in contrast to XWi which takes almost a con- 
stant value in a wide temperature range, < T < 1. 

Figures 1 (c) and 1 (d) show the momentum-resolved 
magnetic susceptibilities xp (k) and xvv (k) evaluated at 
T = 0.1, respectively. The major contribution to xp 
comes from the vicinity of three Fermi surface sheets, 
while xvv comes from a wide area in a Brillouin zone 
(e.g., around (7r,0) and (0, if) ), where either the d yz or 
d zx orbital is occupied, and the other orbital is empty. 
We note that xvv is brought about by the hybridization 
of the d yz and d zx orbitals due to the magnetic field par- 
allel to the z-axis. In particular, sharp peaks are located 
at (pxiPy) ~ (±0.627r, ±0.627r), where these bands cross 
at the Fermi level. 




Fig. 1. (Color online) (a) Temperature dependences of xp ant ^ 
Xvv in non-interacting case and A = 0. (b) Fermi surface for 
A = 0. (c) X p(k) and (d) xvv(k). 



Next, let us consider the effect of electron interaction. 
In Fig. 2, we show the xp and xvv divided by their non- 
interacting values Xp°' ) and Xvv- Here, we plot xp and 
Xvv by varying U, with U'/U and J/U fixed. Figures 
2(a)-2(d) show that both xp and xvv tend to increase 
with U ; however, their growth rates are quite different. 
Although xp is substantially enhanced by electron inter- 
action, Xw/Xw remains ~ 1.1, at most, in the inter- 
action range considered here. Moreover, Figs. 2 clearly 
show that the inter-orbital repulsion U' and Hund cou- 
pling J affect xp and xvv, quite differently. Figures 2(a) 
and 2(b) show that xp monotonically increases with J, 
while it is suppressed with increasing U' . In contrast, 
Figs. 2(c) and 2(d) show that xvv decreases with J, 
while it grows with increasing U' . The contrastive behav- 
iors of xp and xvv can be attributed to the difference in 
the way electron correlation affects spin and orbital fluc- 
tuations. In multi-orbital systems, a magnetic moment 



can be decomposed into the spin part and the orbital 
part, as M z — l z + 2s z . Without spin-orbit coupling, xp 
(xvv) is equal to 2s 2 (l z ) divided by the applied mag- 
netic field. Accordingly, the magnitude of xp (xvv) is 
affected by a spin (orbital) fluctuation. 

Generally, the intra-orbital repulsion U and the Hund 
coupling J stabilize the high-spin states and enhance spin 
susceptibility. Accordingly, xp grows with increasing U 
or J. On the other hand, the inter-orbital repulsion U' 
stabilizes the orbital moment by prohibiting two elec- 
trons occupying different orbitals at the same site. The 
Hund coupling J also destabilizes the orbital moment by 
facilitating the simultaneous occupancy of different or- 
bitals. As a result, xvv grows with increasing U', while 
it decreases with increasing J. 

These contrastive correlation effects naturally lead to 
the different enhancement rates of xvv and xp noted 
above. Although xp is enhanced by a large intra-orbital 
repulsion U, the increase in xvv is mainly brought 
about by a smaller inter-orbital repulsion U'. Accord- 
ingly, xvv /Xvv becomes relatively small compared with 
Xp/Xp^- We plot the ratio Xvv/xp i n Fig. 2(e) for sev- 
eral values of U'/U, under the relation U' = U — 2J. 
Evidently, the correlation enhancement of xvv is smaller 
than that of xp for a wide parameter range. 
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Fig. 2. (Color online) (a) (b) U dependence of Xp/Xp ■ ( c ) (d) 
The U dependence of Xvv/Xyv' ^ e ^ x J/U = 0.0 for (a) and 
(c), and U'/U = 0.6 for (b) and (d). (e) Ratio of xp to xvv for 
the interaction parameters under the relation U' = U — 2 J. We 
fix A = in (a)-(e). 



Next, we will consider the correlation effects on xp 
and xvv under a finite spin-orbit coupling. In Fig. 3(a), 
we plot the A dependence of xvv for several U's, with 
U'/U = 0.4 and J/U = 0.3 fixed. With this choice of 
U' and J, xvv is only slightly enhanced by inducing U 
for A = 0, consistent with Fig. 2. However, Fig. 3(a) 
clearly shows that xvv increases with U for moderate 
A. xvv/xp increases from 4.5 to 6.3 at A = 1.2 while 
sweeping U from 0.0 to 3.0. 
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To elucidate the origin of this marked enhance- 
ment, we introduce the spin fluctuation x S an( l or_ 



bital fluctuation \ L as X S = jf J cIt{S z (t)S z ) and 

o 

P 

X L = jf J ' cLt{L z (t)L z ), respectively, with S a {L a ) = 
o 

J2i, vv '( r )\ s i( l i)\^) c l v c iri'- In particular, X L corresponds 
to the fluctuation between d yz and d zx orbitals, which is 
essential to xvv at A = 0. We plot xvv together with x S 
and x L m Figs. 3(b)-3(d) with varying J/U. As Figs. 3(c) 
and(d) show, the correlation effects on spin and orbital 
fluctuations are not sensitive to A. \ S (x L ) monotoni- 
cally increases (decreases) with J, consistent with the 
view that spin (orbital) fluctuation is enhanced (sup- 
pressed) by Hund coupling. In contrast, xvv shows a 
non-monotonic A dependence. For a small A (A < 0.2), 
Xvv decreases with J/U, consistent with the case of 
A = 0, whereas, for an intermediate A (A > 0.2), xvv 
grows with J/[/, 21 ) as is clearly shown in Fig. 3(e). 

This non-monotonic behavior of xvv can be ascribed 
to the mixing of spin and orbital degrees of freedom 
due to spin-orbit coupling. Namely, the spin moment 
does not commute with the Hamiltonian for A ^= 0; 
hence, the spin degree of freedom also contributes to 
Xvv- We define the spin- (orbital-) dominant region by 
the criterion Xw[J/U = 0.3] - Xw[J/U = 0.0] > 
(< 0), and show the two regions in Fig. 3(a). Figure 3(a) 
clearly shows that xvv is strongly enhanced in the spin- 
dominant region. Namely, in the orbital-dominant region, 
Xvv is enhanced mainly by an inter-orbital electron in- 
teraction, whereas, in the spin-dominant region, a large 
intra-orbital electron interaction contributes to XWi an d 
Xvv becomes strongly enhanced. Our current analysis is 
based on the band structure of Sr2Ru04; however, we 
confirm that our results are generic in ti g systems by ob- 
taining qualitatively the same behavior for other band 
structures. The details of our analysis will be reported 
elsewhere. 

Here, let us discuss our results, in association with ex- 
periments. By adopting the band structure for Sr2Ru04, 
we revealed that xw/xp is ~ 0.5 in a non-interacting 
case. The electron interaction tends to make this ratio 
smaller for A < 0.2. Actually, the spin-orbit coupling of 
Si'2Ru04 is rather weak. Therefore, it is reasonable to 
conclude that xp dominates xvv for Sr2Ru04, i.e., xp 
has a dominant contribution to the Knight-shift signal. 

On the other hand, for most Ir compounds, spin-orbit 
coupling is estimated to be fairly large. In light of our 
analysis, xvv is highly enhanced by electron correlation, 
if the spin-orbit coupling is so large that the spin de- 
gree of freedom contributes to xvv- This gives a possible 
clue to the anomalously large residual susceptibility of 
Ir compounds. We note that some of the Ir compounds 
are considered as Mott insulators, which are outside the 
scope of our current analysis based on a perturbative 
method. Nevertheless, we consider that our mechanism 
is also relevant to the large residual susceptibility of such 
compounds. It is an interesting future study to extend 
our analysis to the vicinity of metal-insulator transition. 

In summary, we have studied the effects of electron 
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Fig. 3. (Color online) (a) A dependences of XVV for U' /U = 0.4 
and J/U = 0.3. Inset: xp - The arrow shows the transition to the 
band insulator, (b)-(d) J/U dependences of XVVi X L and X S a * 
U = 2.0. They are normalized with the values at J = 0. (e) Xvv 
at U = 2.0 for small A. U'/U = 0.4 is fixed for (b)-(e). 

correlation on the Pauli susceptibility xp an d the Van 
Vleck susceptibility xvv on the basis of the multi-orbital 
Hubbard model. We adopt DMFT combined with IPT, 
and calculate xp an d XWj following the definitions in- 
troduced by Kontani and Yamada. As a result, we found 
that the correlation enhancement of xvv is rather small 
for a small A. Meanwhile, a substantial increase is found 
for an intermediate A, where the spin degree of freedom 
contributes to xvv- 
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